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1. INTRODUCTION 


This report presents the results of a laboratory experiment in convection 

flow through a phase change, across which the density changes and latent 
heat is released. The working fluid is a nematic liquid crystal that under- 
goes a first-order phase change at convenient temperature and pressure. The 
addition of a first-order phase change to the convection flow makes the problem 
much more interesting, because in addition to the thermal buoyancy gpAT there 
are two more sources of buoyant acceleration, i.e., the latent heat release 
and density difference across the phase boundary. Thus, in addition to the 
ordinary Rayleigh number R^, * ATfig d 3 /vic, there are two other, very similar 
parameters. When the imposed temperature difference, AT, is replaced by the 

temperature that is characteristic of the latent heat release, L/c , the 

P 

governing parameter is R T = (L/c )Pg d 3 /vic. And when the density difference 

" p 

across the phase boundary, Ap, is used, the parameter is R^ p * (Ap/p) g d 3 /vx. 

In the above, 3 is the expansion coefficient, g is gravity, d is the layer 
depth, v and tc are the diffusivities of momentum and heat, L is the latent 
heat, and c p is the specific heat at constant pressure. 

The phase-change convection problem is further complicated by the fact 
that there are two possible geometries, even in the thermally unstable case 
(lower boundary hotter). When the lighter phase is present near the lower 
boundary, increases of both R^ and R^ p are destabilizing; i.e., an increase 
in either parameter is similar to an increase in R^ and brings the flow closer 
to the point of marginal stability. The latent heat effect is due to the 
release of latent heat in the hot updraft regions, and the density effect 
occurs because the basic state is statically unstable, i.e., the Rayleigh- 
Taylor geometry. When the heavier fluid lies near the hotter lower boundary, 
increasing R^ is now stabilizing, since in the updrafts latent heat is taken 



in. Increasing R^ is destabilising, since the phase boundary sores opposite 
to the vertical motion so that gravity, tending to flatten out the phase 
boundary, induces upward motion in the region of updraft. The latter effect 
vas first shown by Busse and Schubert (1971), who demonstrated that the 
phase boundary in this case vas not the barrier to motion that had been 
previously supposed (Khopoff, 1964). In the preceding descriptions, it is 
assumed that the latent heat is released vhen the fluid goes from the lighter 
(higher temperature or lover pressure) phase to the heavier one, so that the 
heavier-above case applies to laboratory experiments vhere compressibility 
effects are negligible and the heavier-below case to the geophysical flows 
where such effects dominate. 

Convection flows in which a fluid is thermally driven (Rayleigh-Benard) 
and undergoes a phase change or in which two immiscible fluids are driven 
by a density discontinuity (Rayleigh-Taylor) have a number of interesting 
applications to geophysics, meteorology, and astrophysics. Thermal convection 
through a phase boundary has been used to look at the problem of mantle con- 
vection that occur 8 in the presence of an olivine-spinel phase change at 
U00 km (Schubert and Turcotte, 1971). The unstable two- fluid system has 
been used to explain the formation of geological structures such as salt 
domes (Whitehead and Luther, 1975 ). in meteorology the latent heat release 
due to phase change is a crucial part in the description of cumulus convection 
and in cellular convection over the ocean, especially the closed-cell hexagonal 
convection occurring in marine stratocumulus layers (Schubert, et al. , 1979 )• 

In astrophysical flows the phase change process can be used to model convective 
flows in which heat is released frem nuclear reactions or fluid properties 
changed due to ionization (Spiegel, 1972). 
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The experiments reported here are necessarily limited to the heavier- 
above case, so the lighter phase must be present at the hotter lover boundary. 
Therefore, the conditions that can be examined are: thermal convection in 

each phase separately, convection with negligible thermal gradient and unstable 
phase boundary, and thermal convection including phase charge. The goal of 
the work is to analyze the unique laboratory fluid flows obtained in the pre- 
sence of a phase change. We vill also see if the results, or other experiments 
using these fluids, can be used to understand geophysical flows with phase 
change. 

We begin with thermal convection in each phase, also giving a brief 
review of the fluid dynamics of nematics, and then proceed to consider in 
turn each of the other cases possible in the laboratory. 

2. THERMAL CONVECTION IN A SINGLE PHASE 

The working fluid in the experiments described below exhibits a first- 
order phase change, i.e., one in which the density changes and latent heat 
is released. Above the phase-change temperature the material is an ordinary 
isotropic fluid, and below it is a nematic liquid crystal. Since we are 
interested primarily in the phase change, the experimental configuration and 
operating conditions vill be chosen to minimize the effect of the unique 
anisotropies of the nematic phase. The single-phase experiments are pre- 
sented both to provide a baseline with which to compare the two-phase results 
and to illustrate that the nematic anisotropies do not greatly affect the 
observed flows. 

2.1 Nematodynamics 

A nematic liquid crystal is an anisotropic liquid made of long, rod-like 
molecules that can be aligned along one direction, designated by a unit vector n. 
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No change in properties is observed when n -+--n. When aligned, the direction 
of molecular orientation is easily determined, since the material is birefringent. 
The thermal and electrical conductivity, index of refraction, and viscosity all 
depend on the orientation of n relative to the direction of propagation. Unlike 
Isotropic fluids, nematics can resist torques elastically. The amount of static 
torque present in the fluid depends on the curvature of n. The molecules can 
be aligned by imposing an electrical or magnetic field, by suitable treatment 
of the confining boundaries that is coupled to the interior by the static torques, 
or by shears in the fluid. It is the last method of orientation that will be of 
interest in the present case. 

The stress tensor of nematic fluidB contains four viscosity coefficients 
in addition to that of ordinary isotropic fluids. It is also a function of n 
(see deGennes, 197*0. The method of solution is to solve the ordinary Navier- 
Stokes equations for viscous flow, using the more complicated nematic stress 
tensor, together with an additional differential equation for n that is the re- 
sult of adding the torques on the molecules due to imposed fields, static torques, 
and fluid shears. When there are only small deviations in 11 from the value de- 
termined by an imposed field, the analysis is considerably simplified, and a 
number of different flows have been solved. Shear instabilities (Dubois-Violette 
and Manneville, 1978) and thermal convection (Guyon and Pieranski, 197*0 in thin 
horizontal layers with n parallel to the confining boundaries both give results 
that are unique to nematicB. Thermal convection in a thin vertical slot when n 
is normal to the confining boundaries has been solved by Horn, et al. (1976) for 
the case when the orientation is due to static torques and by Fitzjarrald and 
Owen (1979) when an electrical field is imposed by a feedback circuit. The key 
to these successful solutions has been that the fluid exists in a thin layer so 
that there are only small deviations from a constant imposed n. 



$y contrast, the present experiments are in thick layers, i.e., too thick 
for orientation at the boundaries or external fields to affect the interior of 
the fluid. Thus, the only orienting mechanism left is that of shear, and the 
viscous torques (see deGennes, 1974) vill align the molecules so that n is 
parallel to the direction of motion, i.e., normal to the shear. The primary 
effect of this orientation vill be to specify the effective viscosity coefficient. 

In simple shear flow there are three principal viscosities in a nematic, de- 
pending on the angle between n and the shear. These are illustrated in Figure 1. 
Near the phase-change temperature a typical room-temperature nematic v is approx- 
imately the same as v^, the isotropic value, while is about 25£ less than this 
value and v is nearly double. So at the onset of motion in a randomly oriented 
nematic, v is slightly higher than in an isotropic fluid, since it is the average 
of the values for the three different orientations. As the motion starts, fluid 
shears are effective in orienting n, so that the value of viscosity becomes very 
close to V^. Therefore, there must be a slight decrease in viscosity as the 
motions 8 tart. 

The thermal diffusivity is also anisotropic in a nematic, with the value 
of K| - 1.5 for room-temperature nematics. Here tc 1 1 and ic^ are the diffusiv- 
ities along and normal to n, respectively. It is clear that, because of the two 
anisotropic effects — lowering of the effective viscosity and heat focusing by 
the anisotropic thermal diffusivity — the onset of motion in a convective cell 
will be considerably different in a nematic than it is in an ordinary isotropic 
fluid. There should be a hysteresis, with motion continuing at lower than that 
necessary for initiation of motion. Once the motion has started with R^ greater 
than the Initiating value, however, the flow should be very similar to that of the 
isotropic phase. In this case the viscosity becomes very nearly v^, and since the 
convective heat transfer is greater than the conductive part, the conduction 
anisotropy is not a dominant effect. 
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In the experiments to follow, the peculiarities of nematics will be avoided 
hy concentrating on the finite-amplitude flows that exist once the flow is veil 
established. It may be noted that the discussion given previously for a randomly 
oriented nematic is considerably different than for a well-oriented one. 'Then 
the molecules are Initially oriented by an external field or surface treatment 
so that n is parallel to the confining plates, the motion starts at R^-l instead 
of R.-1700 for an isotropic fluid, because of the heat-focusing effect of the 
conductivity anisotropy. Coupling of the heat focusing to the thermal convec- 
tion by means of changes in n can also lead to oscillatory, or overstable con- 
vection because of the difference between the convection and orientation time 
scales (Guyon, et al., 1979). And when n is initially normal to the plates, 
convective motions can occur when the fluid is heated from above (Guyon and 
Pieranski, 197*0 . 

2.2 Description of Experiment 

A sketch of the apparatus is shown in Figure 2. The fluid is confined by 
means of glass plates in a horizontal layer *45 mm diameter by 6 . k mm deep. 

Water is circulated to keep the confining plates at constant temperature. Con- 
ductivity of the glass is approximately 6 times that of the worh ng fluid, and 
the vater flow rates are large enough so that nearly an isothermal boundary 
condition is maintained. The working fluid is a cyanobiphenyl, K-15 from BDH 
Chemicals. It is nematic from 22.5°C to 35°C and is very stable, being resis- 
tant to degradation by moisture, exposure to the atmosphere, and uv radiation. 
Unfortunately, not all the physical constants are known. The values Bhould 
not be significantly different from those of other room-temperature nematics 
(such as MBBA — more widely studied, but much less stable), so that estimates 
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65 mm dia glass plates 



Figure 2. Sketch of convection cell. 






based on such typical values should be applicable to K-15. The physical 

propertiei of interest that are known for K-15 are: 

p • 1.02 gm/cm 3 i 25°C 

0 • lO -3 °C" 1 6 25°C 

n /n - 1.7/1. 5 6 25°C 

e o 

L * 0.8 cal/gm, 

where n e and n o are the extraordinary and idinary indices of refraction. Prop 

erties of MBBA that should be approximately the same as K-15 are: 

2 

. 0.3 cm /s 

„ 10” 3 cm 2 /s 

Ap/p . 2.10' 3 

c - 0.A cal/gm°C. 

P 

We see that Ap/p and L/c^ both are approximately equal to a 2°C change in fluid 

o 

temperature. Using the preceding values, we find that AT * 2 C is enough 
to make 1^.1700 and initiate thermal convection. We find that motion starts 
with AT_1,5°C in the isotropic p’.ase just above the transition, so that the 
above valueB must be close to the true valueu. We intend to avoid the region 
near critical R^, so that a more precise value of the constants is not nec- 
essary. 

2.3 Thermal Convection in One Phase Alone 

With only one phase present and an unstable temperature gradient, the 
famllla'' Rayleigh-Benard results are obtained. Convection rolls are observed 
when AT exceeds some critical value, and the rolls are circular because of the 
circular symmetry of the experiment with downflow near the outer edge. Upflow 
occurs at the center of the experiment with 3 waves across the diameter, corres 
ponding to a wavelength X » 2.3 d. However, the limited abpect ratio and the 
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necessity of having an integer number of waves in the experiment precludes 
making any general conclusions from this measurement of the wavelength. Flow 
visualization in the clear isotropic (lighter) phase is by shadowgraph and 
interferometer that show the temperature distribution in the layer. In the 
liquid crystal (heavier) phase, the birefringence allows observation of the 
fluid motions. As soon as motion starts, the molecules line up along the 
direction of motion due to the action of shear torques on n. ' r hus, the fluid 
appears clear in the regions where it is moving directly toward or away from 
the observer, i.e., the hottest and coldest part of the convection rolls. In 
the region where the fluid is moving ncrr’cJL to the direction of viewing, the 
fluid appears opaque. 

A photograph of the transmitted light is shown in Figure 3 for the condi- 
tions when AT * 2°C and the hotter temperature is Just below the phase-change 
temperature, T c > Upflow at the center and second ring with downflow at the 
first ring and the outer edge is confirmed by direct observation with a stereo 
microscope. With the same AT and the coldest temperature Just above T c , i.e., 
in the isotropic phase, the thermal pattern is seen to be Just the same as for 
the heavier phase. A photograph of the interference pattern is shown in Figure 
U, where the wider spacing between bright lines indicates colder temperatures. 
Careful observation of the interference pattern (and the associated shadowgraph) 
indicates that the diameter of the coldest part of the convection roll is Just 
the same as for the other phase. A roll-like pattern is always observed in 
thermal convection in Just one phase. When the rotational symmetry of the 
flow is disturbed by inverting the cell and shaking it, the roll pattern that 
occurs once the layer is returned to horizontal is no longer circular. Sinuous 
rolls with spacing X_2d are observed that are. similar to the photographs of 
Willis and Somerville (1972) for convection in a rectangular geometry. 
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ri . , Thermal convert io n ~ t he Isotropic phase alone. AT - * . — 

r ‘ " between the interference lines Indicate, cooler temperatures. 
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The circular thermal pattern of the con/ection rolls is first seen in the 
lighter phase when AT£1.5°C. When AT = 2°C, the rolls are vigorous and the thermal 
pattern is quite distinct. In the heavier liquid-crystal phase the roll pattern 
is seen to start as AT is raised above 0.5°C and to persist once Btarted until 
AT is lowered to 0.3°C. It appears, therefore, that the anisotropies of the 
liquid-crystal phase have lowered the critical temperature slightly, but the 
two visualization techniques do not necessarily have the same threshold, so a 
definite conclusion cannot be made. The hysteresis that occurs in the liquid- 
crystal phase is very slight, being almost within the precision of the tempera- 
ture measurement . 

Certainly, the effects of anisotropy are not great and do not result in a 
qualitatively different flow as is the case for uniformly oriented liquid crystals. 
With AT * 2°C the convection flow is the same in either phase, as nearly as can 
be determined. It Is thus inferred that the liquid-crystal anisotropies have 
only a slight effect on the fluid motions to be described below. 

3. RAYLEIGH- TAYLOR INSTABILITY IN THE TWO-PHASE SYSTEM 

A statically unstable condition can be produced in a nearly isothermal 
layer, either by slightly elevating the lower temperature when the whole fluid 
layer is in the heavier phase Just below T c> or by slightly lowering the upper 
temperature when the whole fluid is Just above T c> In the first case a thin 
layer of lighter fluid forms at the bottom and in the second a thin layer of 
heavier fluid forms at the top. As the depth of the thin layer gets larger, 
the interface is seen to become distorted, finally leading to spouts of fluid 
shooting out of the thin layer. The spout of unstable fluid thus formed be- 
comes the center of a cellular pattern that continues to have a convective 
flow. The flow is quasi-steady with new cellB starting up and old ones dying 
out but always with approximately the same spacing. 
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Figure 5 is a photograph of the transmitted light of the cellular pattern 
formed when the lower temperature is 0.1°C above T c - The bright circles are 
regions of upshooting out of the layer cf clear, light fluid near the lower 
bcmdary. The slow return flow is along the boundaries of the cell, A parcel 
<■; fluid thus starts in the light phase, shoots up in the center, spreads out 
along the slightly cooler top boundary changing to the heavy phase, and then re- 
turns to the center of the cell to commence another trip. Figure 6 is the cor- 
r* sponding photograph of the pattern when the upper boundary is 0.1°C lower than 
T , the darker circles now being the downshooting heavier fluid at the center of 
the cells. 

A sketch of the flow patterns is presented in Figure 7, based on observa- 
tions with the stereo microscope. It is seen that sinking motion is associated 
with changing from heavy to light (taking in latent heat), and rising motion 
with changing from light to heavy (giving off latent heat). This behavior is 
as expected in the qualitative discussion in the introduction and is confirmed 
by watching t v . • phase-change process. When AT = 0 and a parcel of fluid is 
changing from heavy to light (melting), it always sinks, and when it is changing 
from light to heavy (freezing), it always rises. Thus, spheres of heavy fluid 
that are becoming larger rise, and those that are becoming smaller sink, when 
AT = ). 

The Jii’tiation of the cellular motions is described by the analysis of 
Whitehead and Luther (1975). who have analyzed the linear and finite-amplitude 
R-yieigh-Taylor problem. The linear results for p thin layer of unstable fluid 
when the viscosities are equal are that the wavenumber of maximum growth is 
k ■ £tt/A = 1.8/2h, and the growth rate is n = 0.15 g (Ap/p) h/v, where h is the 
dep,'j <-f the thin layer. In the present ease the calculated time for growth Is 
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Figure 5. Buoyant motions due to unstable interface. Upshooting plumes from the 
thin layer of light fluid near the lower boundary create cellular 
motion. AT ■ 0. 1°C. 
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Figure 6. Buoyant motions due to unstable interface. Downshooting plumes from 
the thin layer of heavy fluid near the upper boundary create 
cells with downgoing centers. AT = 0.1°C. 
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5 oln when h * 0.1 mm and 30 s when h = 1 mm. Thus, when the layer depth hr," 
grown to 1 mm, the instability should grow quickly, with a horizontal wavelength 
approximately equal to the depth of 6 mm. These are roughly the same values 
observed for both the upshooting and downshooting cells. The finite-amplitude 
analysis of Whitehead and Luther indicates that the planform of the motion 
should be hexagonal, with the sign of the motion so that the fluid Jets out of 
the thin layer. This result is also confirmed by the present results. 

Once the motion has been established by the unstable density interface 
( cor~esponding to the parameter discussed in the introduction), it continues 
without stopping. This is in contrast to the classic Rayleigh-Tavlor problem 
where the fluid spout necks down and becomes a bubble at the other boundary. 

The continuing flow must be due to the effect of latent heat release, which is 
of the proper sense and magnitude. But the fact that this latent-heat driven 
convection needs the finite- amplitude Rayleigh-Taylor flow to kick it off and 
a small temperature gradient to keep it going would surely complicate analysis 
of the flow. 

4. COMBINATION OP RAYLEIGH-TAYLOR AND RAYLEIGH-BENARD INSTABILITIES 

When an already-existing thermal convection flow in a single phase is al- 
tered by slightly changing the mean temperature (keeping AT the same) so that 
the phase change occurs at one of the boundaries, the planform of the convection 

is immediately altered. With the lower plate temperature Just above T , the 

c 

flow changes from rolls to cells with upshooting centers. This condition is 
shown in Figure 8 where AT * 2°C. Parcels of fluid go right through the phase 
boundary as discussed in Section 3» except with much greater speed. Because of 
the circular symmetry of the initial thermal convection flow the upshooting 
parts of the cells are arranged around the circle of up-flov and the center. 




Figure 8. Thermal convection with phase change. Lower boundary is just 
above the phase-change temperature. Cells have upgoing 
centers. AT = 2°C. 
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When the upper plate temperature is Just below T c , the cells have downshooting 
centers and are arranged around the downflow ring of the initial flow. This is 
shown in Figure 9» where the dark line is approximately 2d. In both cases, the 
flows are steady state. 

The changing of the convection planform from rolls to hexagons when some- 
thing is present to disturb the vertical symmetry of the flow is akin to the 
work of Busse (1962) for vertical variations in properties and Krishnamurti 
(1968, 1975) for time variations and boundary suction. In each case the verti- 
cally symmetric rolls must give way to the asymmetric hexagons. When convec- 
tion occurs in a liquid, for example, the hexagons must have upshooting centers 
so that the Jet occurs in the low- viscosity, warm fluid. Gas hexagons go the 
other way, because the low viscosity is in the cold fluid. In the present 
case the Rayleigh-Taylor instability provides the asymmetry to cause the change 
in planform. The latent heat release is available to make the motion go faster, 
but since it is present in both updraft and downdraft, there is no contribution 
to asymmetry. 

5. SUMMARY AND CONCLUSIONS 

The laboratory experiments described here have explored the effect of phase 
change on buoyancy-driven flows in a horizontal layer. Thermal convection in 
each phase separately gave a basis with which to compare the two-phase results 
and showed that the anisotropies of the liquid-crystal working fluid did not 
greatly affect the results. When no appreciable thermal gradient was present 
and the phase boundary was unstable, the experimental results correspond to 
the analysis of Whitehead and Luther (1975). The continuation of flow after 
its Initiation by the unstable phase boundary is a new result and is unique 
to these experiments. The release of latent heat as the fluid passes through 







Figure 9. Thermal cuuvecuun mu pua»c um^c, uuunaary is just below 

the phase-change temperature. Cells have downgoing centers. AT ■ 

2°C. Length of dark bar is approximately twice the fluid depth. 
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the phase boundary must be the primary cause of this quasi-steady flow. 

When thermal convection exists in the presence of a phase change, the plan- 
form was observed to be cellular, instead of the rolls in thermal convection 
alone. The cause of this change in planform is the vertical asymmetry due 
to fluid shooting out from the thin layer across the density discontinuity 
of the interface. 

The results have shown that a steady convection flow can exist right 
through a phase boundary, and that the flow is profoundly changed by 
presence of a change of phase. The limits of the laboratory and comp . JU y 
of the problem limit the applicability of these results to phase-change flows 
in geophysics, meteorology, and astrophysics, however. The geophysical case 
that may be applicable to mantle convection, where the heavy fluid lies below, 
cannot be attained in the laboratory. Geological flowr such as salt domes 
and volcanism are perhaps better modeled by fluids with vastly different viscosi- 
ties, as used by Whitehead and Luther. The meteorological case of cumulus con- 
vection driven by latent heat release is certainly an important one and has been 
the object of previous laboratory experiments. In such a case, the release of 
latent heat occurs only in the updraft regions, which is different from the 
fluid used here. 

In marine Btratocumulus layers that are driven by latent heat release the 
density interface, i.e., the boundary between the marine layer and the free at- 
mosphere, is a stabilizing influence (in contrast to the mantle convection where 
latent heat is stabilizing and the density interface destabilizing). So the 
present results do not apply, because in each of these cases the principal 
driving forces sure in the wrong sense relative to each other. It appears that 
the fluid used in these experiments cannot be used to directly model 
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geophysically relevant experiments. However, the results do show that the latent 
heat released in the phase change can keep the fluid l^j^ic^ion when it would 
ordinarily die out; and the density difference across the phase boundary can 
change the usual convectiou rolls to cells. Both these effects are similar 
to those that likely occur when ,jhase change is present in geophysical or 
meteorological convection flows. 
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